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ON THE APPLICATION OF THE METHOD OF LEAST SQU'S 

TO THE BEDUCT'N OF COMPARISONS OF LINE-MEASURES 

AND TO THE CALIBRATION OF THERMOMETERS. 



BY T. W. WEIGHT, B. A., C. B., DETROIT MICH. 

1. Comparisons of line- measures and observations for the calibration 
of thermometers may be discussed together, as there is no essential difference 
in the mode of making the observations or of reducing them. 

Let AB be the line-measure divided into equal A B 

parts as nearly as maybe at the points 1, 2, 3, ... '„ \ ', 'a '4 
The problem is to find the corrections to the grad- C D 
uation marks at the several points of division. In ' ' 
making the necessary observations, microscopes are mounted over the points 
and 1 and readings made. The bar is then moved along until the points 
1,2; 2,3, .. . are under the microscopes and readings are made each time. 
For the purpose of facilitating the work an extra scale CD may be employed 

and comparisons made between it and each of the spaces 1, 1 2, in 

succession. The correct's then are the amounts by which the marks should 
be changed in order to be in their true positions. When the corrections 
have been applied to the distances 01,12,.., these distances are all reduced 
to the same length. 

With a thermometer AB, if 0, 1, 2, ... . denote graduation marks, then 
by breaking off a column of mercury approximately equal in length to the 
distance 1, we may by reading the ends of this column when in the posi- 
tions 01,12,... determine the degree of uniformity of the bore of the tube, 
whence may be derived at the several points the calibration corrections, as 
they are called. The corrections to the graduation marks for calibration 
when applied to the observed values of the spaces 01,12,... change them 
to what would have been found had the tube been of uniform bore through- 
out its whole length. 



—34- 

2. As the precision with which comparisons of line-measures can now be 
made is very great, it is fitting that the reduct'n of such comparisons should 
be carried out by the method of least-squares. With thermometers however 
such niceness of reading has not yet been attained and approximate methods 
of reduction may consequently be used to advantage. By such methods re- 
sults will be found agreeing as closely with those found by the rigorous 
method as the quality of the observations will warrant. In investigating 
the reliance to be placed on approximate methods we must understand the 
rigorous method to see just what liberties we do take. 

3. In a line-measure AB, let the spaces 1, 1 2, . . . be compared with 
the auxiliary space CD and let Xq, a;j, . . . be the graduation corrections at 
the points 0, 1, . . • Then, for the first comparison, if H^, i?^ be the mi- 
croscope readings at C, D and 31 q, 3f^ be the readings at 0, 1 we have 

where y is the amount by which CD differs from the space intended to be 
represented by CD and v is the residual error of observation. 

Hence, for simplicity taking 4 spaces only, the observation equations may 
be written 

x^—xi —y—{l 2) = «2 { 

a^3— ■'^2 —y—(2S) == vg 

«4— a's —2/— (34) = Vi J 

where (01) = — (Jf^— li'o)-f(i2i— -Bo) 



\, (1) 



As these equations contain 6 unknowns and are themselves only four in 
number, the number of unknowns must be arbitrarily reduced in order to 
carry out the sol tion by the method ot least squares. The simplest sup- 
position is to make a;,, = 0, x^ = 0, y = a, known quantity. All of the 
readings being supposed to have been equally well made and the observa- 
tio equations to be of the same weight, we have the normal equations 

2a; J— x^ = (01)'— (12)' ■) 

_ x^+2x^- x^ = (12)'-(23)' y, (2) 

— x^ + 2xr, = (23)'-(34)' j 

from which the unknowns are found. 

Writing [all, [6?], [cZ] for the constant terms of the normal equations 
and solving 



= i(3 
= i(2 

= i( 



"aZ]-|-2 

"aZ]-f4 

aq+2 



+ W 

H-2[cr 



libq+sici}) 

Hence weight of x^ = f , of x^ = f , of x^ = f . 



(3) 
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4. Generally if the number of spaces is n, 



n.x. 



= (n-l) 

= (71-2) 

= (n-3) 



■«?]+ {n—2)[bl 
'al^ + 2(n—2)lbl' 
'al]+2{n—2)[bl 



+ (n-3) 
+2(n— 3) 
+3(«— 3) 



[c?]+ - 



Hence the weight of any correction as a;„, is 



m{n 



m 



Since this expression remains unchanged when n — m is written for m, and 
is a minimum when m = Jn, it follows that for graduation errors equally 
distant from the ends the weights are equal, and the weight decreases as we 
approach the center where it is a minimum. 

5. It is in accordance with general experience and the principles of least 
squares that instead of spending all of the time of observation on the direct 
measurement of the spaces 1, 1 2, . . . better results would be obtained by 
spending part of it in reading combinations of these spaces as 2, 2 4, ... ; 
3, 1 4, . . . Let then, for uniformity, each space and all possible combi- 
nations of spaces be equally well read. This will require instead of a single 
auxiliary space CD an auxiliary bar divided as nearly as possible the same 
as the bar AB. The observation equations therefore when written in full 
are 



-a; 



*U Q "■ *i/q 



ill 

—2/1 
—2/1 
-2/i 






—X, 



—2/2 
—2/2 



— »A 



—X-, 



—2/3 



-(01) = v^ 
-(1 2) = V, 
-(2 3) = 0,3 
-(3 4) == «4 
-(0 2) = V, 
-(1 3) = V, 
-(24) = V, 
-(0 3) = V, 
-(14) = v. 



-y,-{04) = 



^10 



(4) 



10 equations with 9 unknowns. 

6. Two cases may arise. In the first place the values of the spaces on 
the auxiliary bar CD may be known in terms of the accepted standard. 
We had an example of this on the Lake Survey when the J mm. spaces on 
the Repsold Meter were compared with spaces on the Troughton & Simms 
Inch. The values of the spaces on the inch had been previously deter- 
mined. 

In this case then 2/i) 2/2> S'sj 2/4 ^^^ known and the observation equations 
may be written more simply 
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-a; 3 



-iCo— (0 1) = vi ^ 
-(1 2) = V, 
-(2 3) - ^3 
-(3 4) = v^ 

-x,-{0 2) = «5 
-(1 3) = V, 
-(2 4) = ^, 

-a;o-(0 3) = Vg 



^n 



-(14) 
-(0 4) = 



'10 



(5) 



The normal equations are 



(6) 



ixo—Xi—x^—Xs—Xi = — (0 1)— (0 2)— (0 3)— (0 4) "I 

4x^-x^-x^-x^ = (0 1)-(1 2)-(l 3)-(l 4) I 

4x,-x,-x^ = (0 2)+(l 2)-(2 3)-(2 4) } 

ix^-x^ = (0 3)+(13)+(2 3)-(3 4) | 

40^4= (0 4) + (14)-F(2 4)-f(3 4)J 

Adding these equations there results 

= 0. 

The reason of the indeterminate form is that neither the interval nor the 
terminal point has been fixed. Consequently any arbitrary relation may be 
assumed between the corrections. 

If as above we put a^^ = and x^ = in the observation equations then 
the normal equations become 

4x^—x^—x^ = (01)— (12)— (13)— (14)-) 
4x,~-x, = (0 2)-^(12)-(2 3)-(2 4) I, 
40^3 = (0 3)+(13)-t-(2 3)-(3 4)j 

from which x^, x^, x^ may be found. 

However, by assuming the arbitrary relation 

Xo+x^+x^+x^+Xi = 
between the corrections we obtain a much more simple solution 
adding this relation to each of the normal equations (6) we have 

dx^ = —(01)— (0 2)— (0 3)— (04) = 
5a;i = (0 1)— (1 2)— (1 3)— (1 4) = 
5x^ = (0 2)+(12)— (2 3)— (24) = 
5.T3 = f0 3)-h(13)4-(2 3)-(34) = 
5x^= (0 4)-f(14)+(2 4)-h(3 4) = 

The whole solution may therefore be conveniently arranged in tabular 
form. The sums of the horizontal rows are first formed and then placed 
in the proper vertical columns with the signs changed. The vertical col- 
umns are next added. 



'al'] suppose^ 

dl'] 

"en " 



(7) 



(a) 
For by 



(8) 



—37- 






1 


2 


3 


4 




— Sum^ 


(01) 
— Sum 2 


(0 2) 
(12) 

— Sum 3 


(0 3) 
(13) 
(2 3) 

— •Sum4 


(04) 
(14) 
(24) 
(3 4) 


Sumj 
Sum 2 
Sum 3 
Sum 4 


[aZ] 


[W] 


[cO 


Ldn 


bU 





7. The preceding method of transforming the normal equations (6) into 
(8) by means of the arbitrary relation (a) is allowable. For since the sum 
of the coefficients of the unknowns in each of equations (6) is zero, whatever 
values of x^, x^, . . . satisfy these equations, the values Xf,-^k, Xi-\-^ , . . 
where X is any constant will also satisfy them. Hence whatever set of val- 
ues is taken to satisfy the equations the differences x^ — Xq, x^ — a?o, . . will 
be the same. Therefore by arbitrarily fixing the initial point of gradua- 
tion determinate values are found for the corrections at the other points. 

8. In the second place if, as is usually the case, the values of the spaces 
on the auxiliary bar are unknown the last of equations (4) becomes without 
meaning and must be omitted. This gives 9 observation equations involv- 
ing 8 unknowns from which equations, as the initial graduation point is not 
fixed, the normal equations resulting would be indeterminate. Putting Xq 
= and a;4 = we have 9 equations and 6 unknowns from which follow 
the normal equations 



+ 4a;i— x^— 
— iCj-)-4a;2- 



+ 2/2+ 2/3 = 



aj, Xn -\-4:Xo 



+ 

+ 



42/1 



—2/2— ^3 = 



+32/2 



where 



+2^3 = 
[al] = (01)-(12)-(13)-(14) 



r«^i 


■\ 


[bi] 




[d] 




fdr 


)■ 


fel] 




[/] 


J 



(9) 



This is a more general solution of the problem. 

The elimination of the normal equations from (4) may be considerably 
shortened by an artifice employed by Hansen in his treatise Von der Be- 
stimmung der Theilungsfehler eines gradlinigen Maassstabes, but the amount 
of work is still very considerable. 

9. The Precision of the values of the Unknowns, x^, aij, x^. 

Eliminating 2/2, 3/3 from the normal equations (9) the resulting equations 
may be written 
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(10) 



If the solution is finished by the Gaussian method of elimination the 
weights of the unknowns will be found to be 

weight of x^ = y-, of x^ = ^, of x^ = -*/. 

Generally, if the number of graduation marks is n+1 the weights of the 
corrections to the two outside marks will be found to be n{n-\-l')-r-{2n — 1) 
and of the center marks |-(n+l), where n is an odd number. 

The p. e. of an observ'n of weight unity is found from the usual formula 



.6745 



^■> 



[vv] 



[no. of obs. — no. of indp't unknowns 

where v is the residual error of observation. 

Hence the p. e. of x-^, x^, x^ are known, since their weights are known. 

10. The solution of equations (9) is quite complicated. The amount of 
labor increases very rapidly with the increase of the number of spaces read 
on. Accordingly whenever it is not possible to make use of the simple and 
symmetrical form (8), an approximate form of solution should be employed. 
This is more especially true of thermometers where the carrying out of the 
rigorous solution would entail an expenditure of labor altogether out of pro- 
portion to the increased gain in accuracy. "To devote punctilious attention 
to hundredths of a thermometric degree in results derived from data which 
are uncertain by ten or twenty times that amount is to 'strain at a gnat and 
swallow a camel.'" (Amer. Jour, of Science, Feb., 1882.) 

Instead of finding the corrections to the graduation marks directly, the 
corrections to the several spaces between the graduation marks may first be 
found and thence the cerrections to the marks at the ends of the spaces. 
Thus if Zj, 22> • • • denote the corrections to the spaces, then taking our ex- 
ample of four spaces 

Zi ^^^ X-i x^ I 

}. (12) 



— x^ x^ 



^8 = ^3—^2 j 
2^ = X^ 3*3 J 

Hence by addition 

Z1 + Z2+23+24 = ^4— «o 
= 0. 

From the observation equations (4) we have by subtracting in pairs 
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— H— 


[(1 2)-(0 1)- 


= V 




[(2 3)-(12)" 


= v' 


— 


-(3 4)-(2 3)- 


= d' 


—H— 


-(1 3)-(0 2)- 


= v' 




-(24)-(13)- 


= v'' 


— 2l — 


[(14)-(0 3): 


= V 



= V' '] 



J 

Hence considering (1 2) — (0 1), (2 3) — (1 2), . . . as the independently ob- 
served quantities (which is not correct as they are entangled), and making 
the sum of the squares of v', v'\ ... a minimum we have the normal eq'ns 

] 

^(14) 



32i— Z2— 23— 24=- 


-(12)-(01)] 


- 


L(13H0 2)] 


— 


[(14)-(0 3)^ 


S^a— ^3— 24 = + 


[(12H01)- 


— 


:(2 3)-(12)- 


— 


(2 4)-(13) 


323— 24 = 


;(2 3)-(12); 


— 


_(3 4)-(2 3): 


— 


(13)-(0 2)- 


324 = 


[(3 4)-(2 3)" 


+ 


[(2 4)-(13;] 


+ 


[(14)-(0 3)] 



which when added give = identically. 



But 



+ 22 + «3 + 24 = 0. 



Hence adding this relation to each of the normal equations we have the 
values of the unknowns. Thus 

42, = _[(12)-(0 !)]-[(! 3)-(02)]-[(14)-(03)] 

- - - _ - - - _ (15) 

This form of approximation when applied to thermometers is known as 
Neumann's method. It is quite analogous to the system used by the Prus- 
sian Landestriangulation in the adjustment of horizontal angles and the 
principle involved is therefore really due to Gauss and Hansen. 

For convenience the computation of the corrections may be arranged in 
tabular form as follows : 



2l 


^2 


^3 


24 




— Sum J 


+(12)~(01) 
— Sum 2 


+(13)-(02) 
+(23)-(12) 

— SutUg 


+(14)-(03) 
+ (24)-(13) 
+ (34)-(23) 


Sum, 
Sum, 
Sum 3 


42, 


42^ 


423 424 





The corrections x-, , x. 



are then found from the relations 



= 2- 



^2 — ^i~r^2 
x^ := a;2-r23 

From equations (15) and (16) 

4a?i == —(12 — 01)— (13 — 02) — (14 — 03) 



4x2 = —(13 — 02)— (14 — 03)- 
4x<i = —(14 — 03)— (24 — 1 3)- 



■(23 — 12) — (24- 
-(34 — 23) 



(16) 



■13)Wl7) 
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11. The Precision. — The p. e. of the unit of weight is given by 



r = .6745 



V 



I'o'v''] 



K«— 1)(«-2)' 



(18) 



since the number of equations = Jn(n — 1) and the number of independent 
unknowns = n — 1, the number of spaces being n. 

Also since z^, z^, z^ result from the normal equations with the same w't 
they have the same p. e. From equations (17) it is evident that the p. e. of 
ajj and rcg is the same as that of z-^ but the p. e. of x.2 is i/f as great. 

Generally when the number of spaces is an even number n, the p. e. of a;^ 
is to the p. e. of the center graduation mark as i/(« — 1) : Jn. 

Example. — The following were the observed values of the lengths of the 
45°, 90° and 135° columns of Thermometer Green 4470, made to determine 
the calibration corrections at the 77°, 122° and 167° points. (Amer. Jour. 
Science, Vol. XXI, p. 374.) 



45° col. 


90° col. 


135° col. 


44.68 
.71 

.72 
.75 


90.07 
.09 
.11 


1 
134.61 
.68 



The corrections at the points 32° and 212° are assumed to be zero. 

First Solution. By the method of Least Squares. The observation eq's 
are with the preceding notation, x^, x^, x^ denoting the corrections at the 
77°, 122° and 167° points respectively, 



'^i—yi —0.32 = 





^•2—^1—2/1 —0.29 = 





«3— ^2 — 2/i —0.28 = 





—x^ — )/i —0.25 = 





+«2 —2/2 +0-07 = 





^3 —^1 —2/2 +0-09 = 





~^2 —2/2 +0.11 = 





^3 —2/3—0.39 = 





-•^1 -2/.S-0.32 = 





The Normal equations are 




4x^— X2— a;3+ 2/2+ 2/3 = 


0.33 


— x^+ix^— x^ = 


0.05 


— Xi~ a;2+4a;3— y^— y^ = 


—0,20 


42/1 = 


—1.14 


% — «3 + %2 = 


0.27 


ajg — .Tg +22/3 = 


—0.71 


Solving we find 





-41— 



X, = +0.031 



^2 



2/1 = —0.285 
+0.028 2/2 = 0-090 

+0.031 2/3 = —0.355 

Substitute these values in the observation equations and the residuals re- 
sult. Squaring and adding, 

[««] = 0.0003. 
Hence for the p. e. of a single observation 

r = .6745 x/-^«~|-^ = ±0.0007. 
. • . p. e. of a?! = .007v/.35 =0°.004 






= .007i/.3 =0°.004 



" x^ = .007,7.35 =0°.004 
Second Solution. By Neumann's formula — Arranging according to table 



2l 


^2 


^3 


«4 


Sums 


+.12 


—.03 
+.03 


—.02 
—.01 

+.03 


—.07 
—.02 
—.03 


—.12 
—.03 
—.03 


+ .12 


00 


00 


—.12 





Also 



= +0°.03 
= +0°.00 



= +0°.00 
= — 0°,03 



^1 

x„ 



2l 

Xi +22 

X. 



= 0'>.03 
= .03 
= .03 



r 
and p. 



>3 ^2"T^3 

= .6745v^-^^^ = ±0.007 

e. of .171 = .007i/i = ±0°.004 
■ ' = .007i/i = ±0 .004 

±0 ,004. 






iCg = .007i/J 



THE MULTISECTION OF ANGLES. 



BY PEOF. J. W. NICHOLSON, A. M., LA. STATE UNIV., BATON EOUGE, LA. 

To divide an angle into any number of equal parts, we use an instrument 
called the muUisedor, the simplest form of which is the square sot; which 
consists of two perpendic- p ^ 

ular arms os and ct, the 



length [ct) of the latter be- 
ing equal to the width {cf) 



